The types of Mott insulator 
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There are two classes of Mott insulators in nature, distinguished by their responses to weak doping. 
With increasing chemical potential, Type I Mott insulators undergo a first order phase transition 
from the undoped to the doped phase. In the presence of long-range Coulomb interactions, this 
leads to an inhomogeneous state exhibiting "micro-phase separation." In contrast, in Type II Mott 
insulators charges go in continuously above a critical chemical potential. We show that if the 
insulating state has a broken symmetry, this increases the likelihood that it will be Type I. There 
exists a close analogy between these two types of Mott insulators and the familiar Type I and Type 
II superconductors. 



The study of Mott insulators and doped Mott insula- 
tors constitutes one of the central focuses of condensed 
matter physics of the past two decades. At a simple level, 
a Mott insulator is a system in which the strong repul- 
sion between particles impedes their flow; the simplest 
cartoon of this is a system with a classical ground-state 
in which there is one particle on each site of a crystalline 
lattice and such a large repulsion between two particles 
on the same site that fluctuations involving the motion 
of a particle from one site to the next are suppressed. 
The Mott insulator is thus essentially classical in nature 
(and so accessible to particles of any[l] statistics.) In 
contrast a band insulator [2], including all familiar semi- 
conductors, is a state produced by subtle quantum inter- 
ference effects which arise from the fact that electrons are 
fermions. Nevertheless, one generally considers band in- 
sulators to be "simple" because the band theory of solids 
successfully accounts for their properties; although the 
physical origin of the insulating character of a Mott in- 
sulator is understandable to any child, other properties, 
especially the response to doping, are still only at best 
partially understood. 

There is even debate concerning the definition of a 
Mott insulator. In many of the most familiar cases [3] 
involving electrons, one can, in theory, continuously vary 
the strength of the repulsion between electrons, such that 
in the limit of vanishing strength interactions the system 
can be well understood on the basis of band theory, while 
in the strong interaction limit Mott physics better ac- 
counts for the insulating properties; if there is no phase 
transition between these two limits, it is not clear at what 
interaction strength one should start thinking of the sys- 
tem as a Mott insulator. However, recently an interest- 
ing Bose Mott insulator has been produced by trapping 
bosonic neutral atoms in an optical lattice [4]; there is no 
non-interacting limit in which bosons are insulating, so 
no such ambiguity arises here. 

Mott states, in addition to being insulating, can be 
characterized by the presence or absence of a spon- 
taneously broken symmetry (e.g. spin antiferromag- 
netism), by the nature of the low energy excitation 
spectrum (e.g. gapped or gapless)[5], and by the 
presence or absence of topological order and charge 



fractionalization[6, 7, 8, 9]. To this list, we add a "type" 
index that classifies Mott insulators into two categories 
depending on their response to doping. 

Generally speaking, states with charge gaps (this in- 
cludes both band and Mott insulators) occur in crys- 
talline systems at isolated rational "occupation num- 
bers," v — v* , where v is the number of particles per 
chemical unit cell. By "doping" we mean a process which 
causes the occupation number to shift away from v*. 
When lattice translation symmetry is not spontaneously 
broken, v* is typically an integer for bosons, consistent 
with the simple classical cartoon described above, and an 
even integer for spin 1/2 fermions (e.g. electrons). The 
significance of this latter statement is that a necessary 
condition for an electronic system to be adiabatically con- 
nected to a weakly interacting band-insulator is that v* 
be an even integer. Insulating states can also occur when 
v* is a fraction (for fermions this includes odd-integer 
v*). Usually when that happens, translational symme- 
try is spontaneously broken so that the new unit cell has 
integral occupation number (or even integral occupation 
number for fermions)[10, 11]. For instance, electronic 
Mott insulators with v* = 1 often exhibit antiferromag- 
netic Neel long-range order, which doubles the unit cell 
leading to an effective v*f * = 2. Nevertheless, we will 
mention [6, 7, 8] model bosonic systems for which the 
Mott state can be shown to have no broken symmetries 
for v* = 1/2, and fermionic systems in which v* = 1. 
(Currently, no laboratory system has been found which 
unambiguously exhibits this exotic behavior.) 

The purpose of this paper is to address the response 
of Mott insulators to light doping (i.e. v — > v* — 5). Our 
central observation is that there are two types of Mott in- 
sulators (henceforth referred to as Type I and Type II) . In 
a Type I Mott insulator an increasing chemical potential 
induces a first order phase transition from an undoped 
state to a charge-rich state so that changes occur discon- 
tinuously. In other words, there is a range of "forbidden 
charge density," and two-phase coexistence. In a Type II 
Mott insulator, charges go in continuously above a criti- 
cal chemical potential. Since doping of a band insulator 
is continuous, a doped Type II Mott insulator may be 
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adiabatically connected to a band insulator, but a doped 
Type I Mott insulator is a thermodynamically distinct 
state of matter. 

In the presence of long-range Coulomb forces macro- 
scopic charge inhomogcncity is impossible. However, so 
long as the long-range forces are not too strong, there can 
exist a range of low doping in which the doped charges 
tend to cluster [12]. We still refer to this kind of system 
as a Type I Mott insulator [13]. For instance, if a low 
density of doped charges form puddles or stripes[14], we 
consider this a form of charge clustering, while a Wigncr 
crystal of doped charges is deemed homogeneous. 

What is the situation in real materials? That simple 
semiconductors can be continuously doped, and hence 
are Type II, is well known, although it is not clear they 
are profitably thought of as Mott insulators, even when 
correlation effects significantly renormalize the gap mag- 
nitude. There are many materials in which the insulating 
state is also antifcrromagnetically ordered. It has been 
argued[15, 16, 17] that doping of this class of Mott insu- 
lators generically[18] leads to phase separation, i.e. that 
they are Type I. We will show below that the presence of 
a broken symmetry always increases the tendency of an 
insulator to be Type I, although the conjectured strict 
connection has neither been proven[19] nor disproved. 

For the cuprates and manganates, both of which 
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T = properties of 2D Bose Mott insulators 
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Chemical potential p 
Induced particle density p 
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I. DOPING A TRANSLATIONALLY 
INVARIANT BOSE MOTT INSULATOR 

To begin with, we avoid the interesting complications 
which arise from spontaneous symmetry breaking or from 
the fermionic character of electrons by focusing on the 
simplest kind of Mott insulators formed by spin zero 
point bosons on a lattice at v = 1. This is not a purely 
academic exercise [31], as it applies rather directly to the 



have antiferromagnetic order when undoped, there ex- 
ists considerable evidence that doping induces spa- 
tial inhomogeneity[20, 21, 22, 23, 24]. Conversely, 
Sri_.2;La x Ti03, which is also an antiferromagnet when 
undoped [25], is generally believed [26] to be Type II on 
the basis of the observation of conventional Fermi liquid 
behavior[27] for doping as low as x = 5%, with an effec- 
tive mass that shows a tendency to diverge with decreas- 
ing x; taken at face value[28], this evidence of a quantum 
critical precursor to a continuous metal-insulator transi- 
tion indeed suggests that this material is Type II. 

Another purpose of this paper is to elucidate the 
analogy between Type I and Type II Mott insulators 
and the familiar types of superconductors. In partic- 
ular, when the constituent particles are bosons there 
exists a mathematically precise mapping, the so-called 
"duality transformation" [10, 29], that relates the zero- 
temperature properties of the doped Mott insulator in 
two spatial dimensions, D=2, to the finite-temperature 
response of a 3D superconductor to a magnetic field. Ta- 
ble I summarizes the correspondence between the two. 
This same mapping was exploited previously by Balents, 
Fisher, and Nayak[30] in the context of a theory of a 
proposed "nodal liquid" Mott insulating state; in this re- 
gards, the present paper simply underlines the generality 
of the analogy. 



I 

experiments in Ref. [4]. Consider the following Hamilto- 
nian 



<ij> ~ i 



T > properties of 3D Superconductors 
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where i, j label the lattice sites on a D-dimensional hy- 
percubic lattice, and af creates a boson at site i. The 
first term of Eq. (1) describes the quantum mechanical 
"hopping" of bosons from a site i to its nearest neigh- 
bors j, the second and third terms describe the pair- wise 
interactions between bosons. The U term is a contact in- 
teraction and the Vij terms describe interaction between 
bosons separated by | — r j | . This model is known as the 
extended Bose-Hubbard model, and has been extensively 
studied from various perspectives [32]. 

Let us focus on the U/t — > oo limit. In this "hardcore" 
limit, doubly occupied sites cost too much energy, hence 
are excluded. In this case for large positive [i there is an 
unique ground state 

|Mott>=JJa+|0>, (2) 

3 

in which each site is occupied by one and only one boson 
and hence v = 1. Since this state is separated from all 
other states with the same particle number by an energy 
gap of order U, clearly we have an insulator. (For fi large 
and negative, the groundstate is the empty lattice, |0 >.) 

The behavior of the system upon varying \i (i.e. dop- 
ing) depends on the values of Vij . For Vij = the doped 
holes (i.e the empty sites) interact only through the hard- 
core exclusion. Since the kinetic energy favors a uni- 
form density state with delocalized holes, this limit cor- 
responds to a Type II Mott insulator. Indeed, a small 
concentration of doped holes, 5=1 — v «1, is equiva- 
lent to a system of dilute bosons (defined relative to the 
vacuum state \Mott >) with short-range interactions, a 
problem whose solution has long been known[33, 34]. The 
energy is clearly minimized by delocalizing the holes so 
the ground state in D = 2 is a uniform superfluid with 
superfluid stiffness proportional to tS. Conversely, if Vij 
is attractive (negative), it is clear that the holes will clus- 
ter when | J2i Vij I ^ Thus depending the nature of Vij 
Eq. (1) can describe either a Type I or Type II Mott in- 
sulator. Note that in this example, the Mott state does 
not break any symmetry, so any inhomogeneity of the 
the doped state cannot be attributed to order parameter 
competition. The Type I behavior discussed here is more 
general [35] than that which occurs in the bicritical point 
scenario in the Landau theory of two competing order 
parameters [36]. 

To be still more explicit, consider the case in which 
Vij = —V for (ij) nearest-neighbor sites, and Vij = 
otherwise. Technically in this limit Eq. (1) is equivalent 
to the S=l/2 ferromagnetic XXZ model in a z-direction 
magnetic field 

H = -J xy E {SfSf+SVSp-J, ]T SfSj-h^S?. 

<ij> <ij> i 

(3) 

The mapping between these two models relates J xy to t, 
J z to V, h z to jU + Vc/2 (c = the coordination number). 
The z component of the magnetization is related to the 
boson density according to M z = (v — 1/2)N. (Here 



TV is the total number of lattice sites.) For J z < J xy 
(i.e.V < i), this model has XY order in the absence of 
h z . In this range of parameters varying h z causes the 
magnetization M z to vary continuously, but the ground 
state is uniform and, so long as \M Z \ < N/2, has a non- 
zero component of the magnetization which lies in the 
XY plane and spontaneously breaks the XY symmetry. 
In terms of bosons this means that the insulating state 
is Type II, and the doped system is a uniform superfluid 
for all 1 > v > 0. 

Conversely, for J z > J xy (i.e. V > t) the model is 
effectively an Ising ferromagnet with a fully polarized 
ground state. In this case M z exhibits a discontinuity 
AM Z — N as h z is varied through zero. If we consider 
a constrained groundstate, with \M Z \ < N/2, the state 
exhibits phase separation into two oppositely polarized 
domains. In terms of the bosons, these two domains are 
Mott insulating (v = 1) and empty (v — 0) respectively. 
On the boundary between the two behaviors (J z — J xy ), 
the system is equivalent to a Heisenberg ferromagnet in 
a magnetic field. In the absence of h z the ground state is 
N+l fold degenerate correspond to M z = -N/2, ...,N/2. 
In the presence of h z the state with M z = ±N /2 (± de- 
pends on the sign of h z ) has the lowest energy. For this 
system there is again a first order transition at h z = 0. 
However, precisely at the transition point the boson in- 
verse compressibility k _1 diverges, and a uniform ground 
state exists for any v; at this critical point, only, the sys- 
tem is neither truly Type I nor Type II. 



II. DUALITY TRANSFORMATION FOR THE 
BOSE MOTT INSULATOR IN D=2 

There is a particularly convenient way of thinking 
about Bose Mott insulators that is specific to I? = 2. 
It turns out that for the class of model given by Eq. (1), 
there is an exact mapping, the "duality" transformation, 
that provides us an alternative view of the physics of 
Eq. (1) in terms of the vortices of the boson field[10, 29]. 
It is this mapping that enables us to establish a precise 
connection between the two types of Mott insulators with 
Type I and Type II superconductors[10, 29, 30]. In the 
following we discuss the physical content of the duality 
transformation without going into its technical details of 
the transformation[10]. 

A vortex is a topological defect in the Bose field. When 
a boson is adiabatically transported around a vortex, 
the boson wavefunction acquires a phase factor - the 
Aharonov-Bohm phase 9 — 2ir. In the dual picture, the 
role of the vortices and the particles are interchanged: 
the vortices are the dual particles (they turn out to have 
Bose statistics as well), and when a dual particle is trans- 
ported around an original boson, the wave function ac- 
quires the same Aharonov-Bohm phase[10, 29] discussed 
above. The fact that a boson and a vortex acquire a phase 
when they go around one another implies that bosons 
and vortices can not Bose condense simultaneously. As 
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is well known, in the Bose superfluid phase the vortices 
must be absent or localized, i.e. form an Abrikosov lat- 
tice. Conversely, in the dual phase, where the vortices 
form a superfluid [37], the boson density (and hence the 
dual magnetic flux) must be frozen; the vortex super- 
fluid phase is the Mott insulating phase of the original 
bosons[10, 29]. 

It is important to note that the absence of dynami- 
cal boson density fluctuation is a necessary but not a 
sufficient condition for vortex condensation. A static bo- 
son density acts like a background magnetic field to the 
vortices, which can still frustrate vortex condensation. 
However when the static boson density corresponds to 
an integral u, the vortices see a background magnetic 
flux corresponding to integral number of flux quanta per 
plaquette. (The vortices live on the dual lattice, i.e., the 
centers of the square plaquettes.) This type of flux is 
"invisible" because it can be "gauged away" . The v = 1 
Bose Mott insulator, discussed in the previous section, 
corresponds to precisely this situation. 

When the boson density is a fraction (y = p/q) it 
is also possible for the vortices to condense. This is 
well known [38] in the context of the classical frustrated 
XY model (the XY order parameter is the Bose ampli- 
tude of the vortices) in D — 2. Quantum analogues of 
this, especially for low order rationals (q a small inte- 
ger) can readily be imagined. Such a state is most nat- 
urally accompanied by spontaneous translation symme- 
try breaking[10, 11], as discussed above (and as occurs 
in the classical model), leading to an enlarged unit cell 
with an effective integer v. However, it is possible to 
imagine a more exotic Mott state at v = p/q in which 
the translation symmetry is unbroken. As has been dis- 
cussed in Ref. [6], this could happen if q elementary 
vortices form a bound-state, and these composites then 
condense. Such a composite condensation is unfrustrated 
by an uniform static boson density with v = p/q. More- 
over, since the condensate vortex consists of q elementary 
vortices, charge 1/q bosonic soliton excitation (viewed by 
the condensate vortices as a flux quantum) can become a 
finite energy excitation. Thus there is fractional charge 
solitons! 

In short, a vortex condensate requires the bosons to 
Mott insulate. Consequently we can view a boson Mott 
insulator as a vortex superconductor[10, 29]. Doping 
changes the average background boson density. To the 
vortices this appears as a change in the background mag- 
netic field. The zero temperature boson Mott insulator is 
mapped onto a zero-temperature vortex superconductor 
(with quantum fluctuating two-dimensional electromag- 
netic fields) [39]. However, the quantum partition func- 
tion of a (particle-hole symmetric) [37] two dimensional 
superconductor with a fluctuating gauge field is equiva- 
lent to the classical (i.e. thermal) partition function of 
a three-dimensional superconductor with thermally fluc- 
tuating magnetic field [40]. This leads to the final corre- 
spondence between the Mott insulator and the classical 
fluctuating 3D superconductor summarized in Table I. 



The properties of classical 3D superconductors in a 
magnetic field are generally well known, so most of the 
comparisons in Table I are self-evident. We therefore 
confine ourselves to commenting on a few of the subtleties 
of this comparison: 

Although a Type I supercondunductor typically expels 
any applied field of magnitude H < H c , by choosing 
an appropriate experimental geometry it is possible to 
study its behavior at fixed average magnetic induction - 
which is analogous to studying the Mott insulator at fixed 
doping concentration. This can be done, for example, by 
subjecting a flat slab of a Type I superconductor to a 
perpendicular magnetic field of strength H < H c . It is 
known that when that is done the "intermediate state" 
consists of a mixture in particular, among many possible 
inhomogeneous structures which have been observed[41], 
one of the most common is a laminar (or stripe) structure. 

Much recent attention has been focussed on fluctuation 
effects in the mixed state of a Type II superconductor. 
Mean-field theory predicts that the magnetic induction 
B increases from zero continuously as H is raised above 
the lower critical field H c \. Moreover, the magnetic in- 
duction first appears in the form of flux tubes each en- 
closing a single quantum of magnetic flux which (in the 
absence of disorder) form a regular (Abrikosov) lattice. 
However, thermal meandering of the flux tubes can af- 
fect the physics dramatically near H c \ as well as near the 
mean- field H C 2- Near H c \, the distance between neigh- 
boring flux tubes is much greater than the range of their 
interaction (the London penetration depth A). As a re- 
sult thermal meandering of the flux tubes will melt the 
flux lattice to form a flux liquid [42]. At larger magnetic 
fields the density of flux tubes becomes higher so that 
their interaction can stabilize the flux lattice. This flux 
lattice persists until H — > H c2 where the thermal me- 
andering melt the flux lattice again. The corresponding 
behavior of a Type II Bose Mott insulator is a reentrant 
series of transitions to a superfluid, an insulating crystal 
of doped holes, and again a superfluid as a function of 
increasing doping concentration. 

Finally, it is worth recalling that what determines 
whether a superconductor is Type I or Type II is the ratio 
between the London penetration depth and the core size 
of the vortices, or more physically the sign of the interface 
energy between normal and superconducting regions [41]. 
Type I superconductors have a positive interface energy 
while in Type II superconductors it is negative. As a re- 
sult, the flux tubes effectively attract each other in Type 
I superconductors and repel each other in Type II. As we 
have seen in Eq. (1), this is exactly how we turn a Type 
I Mott insulator into a Type II. 



III. DOPING A MOTT INSULATOR WITH AN 
ORDER PARAMETER 

As we discussed at the beginning of this paper a 
Mott insulating state can be accompanied by transla- 
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tion (and/or other) symmetry breaking. For example let 
us consider Eq. (1) with U/t — > oo and Vij = +V for 
nearest neighbor < ij > and otherwise. For sufficiently 
strong V and fi — the ground state breaks translation 
symmetry and bosons form a checkerboard lattice and 
Mott insulate. In this two-fold degenerate ground state 
the unit cell is doubled. Is this Mott insulator Type I or 
Type II? 

With the above specific choice of U and Vij Eq. (1) is 
equivalent to Eq. (3) with J xy = —t, J z = V and h — 
fi — Vc/2. For this choice of parameters it is known that 
as a function of h z Eq. (3) exhibits a "spin flop" transition 
from the antifcrromagnetic Ising (S z ) ordered phase into 
the ferromagnetic XY ordered phase. Translate this into 
the boson language it implies that it is Type I. 

There is, in fact, a general reason to expect broken 
symmetry to increase the tendency of a Mott insulator 
to be Type I. If there is a broken symmetry in the in- 
sulating state, then there generally is a corresponding 
local order parameter (ip) (for the above example ip is 
the two-sublattice density wave order parameter). We 
now consider the effect of this order parameter on the 
interface energy between the insulating and the metal- 
lic (doped) state. Typically, ip will have a reduced value 
(perhaps 0) in the doped state. Thus, at the interface 
between doped and undoped region a spatially varying tp 
is necessary. The spatial gradient energy 

J rf^lWl 2 (4) 

will thus make a positive contribution to the interface 
energy, increasing the tendency of the Mott insulator to 
be Type I. 



IV. AN EXAMPLE OF TRANSLATION ALLY 
INVARIANT BOSE MOTT INSULATOR AT 

v = 1/2 

A recently discovered system which exhibits a Bose 
Mott insulating state at a fractional v is the quantum 
hardcore dimer model on a triangular lattice. The model 
is defined in terms of bosons (called dimers) which live 
on the nearest neighbor bonds of a lattice, i.e. b\j creates 
a boson on the bond connecting sites i and j. We impose 
a hardcore constraint - this can be thought of as deriving 
from the U — ► oo limit of a contact interaction, which 
forbids two dimers to occupy the bonds emanating from 
a single site, i.e. b\^bijb'- k bjk = for any triplet of nearest 
neighbor sites, < ijk >. Thus, since each bond touches 
exactly two sites, it is possible on any regular lattice with 
one site per unit cell to satisfy this constraint if and only 
if v < 1/2. The simplest (shortest range) Hamiltonian[8] 
one can construct for these dimers is 

H= -tE {im HAikib jk +h.c.} (5) 

+VE {ijk iMAi b ki b ij + ^Ajhjba] 



where (ijkl) denotes a "square" of nearest-neighbor sites, 
such that < ij >, < jk >, < kl >, and < li > are all 
pairs of nearest-neighbors. 

On a square lattice, all zero temperature Mott insu- 
lating phases with v = 1/2 derived from this Hamil- 
tonian arc believed to break translational symmetry, at 
least doubling the unit cell size, so that the effective v 
is integral]?]. However, on a triangular lattice, it has 
now been established [7] that for a range of parameters 
near t ~ V, there is a dimer liquid phase, with exponen- 
tially falling correlation functions and no broken symme- 
tries. From the dual viewpoint, this phase can be viewed 
as a superfluid of vortex pairs [6]. As mentioned previ- 
ously, there are observable topological consequences of 
the vortex pairing, including certain predictable ground- 
state degeneracies on closed surfaces. 

This same model can serve as a paradigmatic example 
of a Mott insulating state of spinfull fermions at an odd 
integer v = 1. To see this, we merely need to define 
fermion creation operators, c\ ^ , for an electron with spin 
polarization a on bond < ij >. In terms of this, we can 
express the dimer creation operators as 



and an additional constraint, ^ CT c\ irj = or 2, 
which can be thought of as a consequence of a strong 
attractive interaction between two electrons on the same 
bond. 

Is this sort of fractional Mott insulator Type I or 
Type II? The intensive study of what happens when 
one dopes such a "spin-liquid" state was initiated by 
the proposal[43] that doping such a liquid would lead 
inevitably to high temperature superconductivity The 
basis for this proposal was the thought that pairing cor- 
relations {e.g. a spingap[9]) might be present already 
in the insulating state, and these would evolve smoothly 
into superconducting pairing upon doping. Central to 
this proposal is the assumption that the spin-liquid is a 
Type II Mott insulator. Indeed, an early triumph of this 
idea was the observation by Ioffe and Larkin[8, 44] that, 
for at least a range of parameters, doped holes in the 
quantum dimer model do not phase separate. Moreover, 
reversing the logic of the previous section, the fact that 
this state has no broken symmetry increases the likeli- 
hood that it is Type II. 

V. CONCLUSIONS 

To summarize, in this work we introduce the concept 
that there are two types of Mott insulator. In particular, 
there exist a whole class of Mott insulators, the type 
I Mott insulators, that become inhomogeneous upon 
doping. The type I behavior does not have to originate 
from order parameter competition. We give a familiar 
lattice boson example where the Mott state exhibits 
no symmetry breaking, however depending on the sign 
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of certain interaction, it can be either type I or type 
II after doping. Through duality these two types of 
insulating states have a close tie to Type I and Type II 
superconducting states. [30] Finally after addressing the 
"zeroth order" physics of doped Mott insulator, we point 
out that many important "first order" issues remain 
unanswered. Among them are "what determines the 
structure of the inhomogeneity in a doped Type I Mott 



insulator?" "What is the ultimate low energy behavior 
of these inhomogcncities?" 



Acknowledgements: We wish to thank C.Nayak for 
useful discussions. This work was supported in part by 
NSF grants # DMR 99-71503 (DHL) and DMR - 0110329 
(SAK). 



[1] In two space dimensions hardcore particles with quan- 
tum statistics between that of fermions and bosons are 
allowed. See, "Fractional Statistics and Anyon Super- 
conductivity" edited by Frank Wilczek, World Scientific, 
1990. 

[2] A band insulator insulates because of Pauli exclusion 
principle. In the non-interacting prototype of this type 
of insulator all available single-particle orbitals are oc- 
cupied. This type of insulator usually possesses a even 
number of electrons per crystalline unit cell. 

[3] For instance, a two leg Hubbard ladder with suitable 
choice of the rung hopping evolves adiabatically from a 
band insulator, for small U, to a so-called Haldane gap 
system, one of the canonical examples of a Mott insulat- 
ing state of electrons without a broken symmetry, in the 
large U limit. 

[4] M. Greiner et al, Nature, 415, 39 (2002). 

[5] A Mott insulator necessarily has a charge gap, however 
it is possible for it to have gapless neutral excitations. 

[6] T. Senthil and M.P.A. Fisher, Phys. Rev. B 63, 134521 
(2001); C. Nayak and K.Shtengel, Phys. Rev. B 64, 
064422 (2001); X-G Wen, cond-mat/0107071; L. Balents, 
M. P. A. Fisher, S. M. Girvin, cond-mat/01 10005. 

[7] R. Moessner and S.L. Sondhi, Phys. Rev. Lett. 86, 1881 
(2001); R.Moessner, S.L. Sondhi, and E.Fradkin, Phys. 
Rev. B 65, 24504 (2002). 

[8] D.S. Rokhsar and S.A. Kivelson, Phys. Rev. Lett. 61, 
2376-9 (1988). 

[9] S.A.Kivelson, D.S. Rokhsar, and J.P.Sethna, Phys. Rev. 
B 35, 8865 (1987); S.A.Kivelson, Phys. Rev. B 39, 259- 
264 (1989). 

[10] D-H Lee and R. Shankar, Phys. Rev. Lett., 65, 1490 
(1990). 

[11] M. Oshikawa, Phys. Rev. Lett. 84, 1535 (2000). 

[12] U. Low et al, Phys. Rev. Lett. 72, 1918 (1994); 

V.J.Emery and S.A.Kivelson, Physica C 209, 597-621 

(1993). 

[13] Strictly speaking, in the limit 8 = v* — v — > 0, since 
the Madulung energy per hole diverges, the groundstate 
in the presence of long-range Coulomb interactions is al- 
ways a Wigner crystal of holes, independently of whether 
in the absence of such interactions the state would be 
Type I or Type II. However, where the Coulomb interac- 
tions are weak, due to the presence of a large dielectric 
constant, clustering occurs in a Type I Mott insulator 
for S > S c , where 8 C is a parametrically small (in powers 
of e~ l ) critical doping concentration. In the absence of 
long-range forces and for small enough 8 , the charges in 
a Type II Mott insulator always form a uniform fluid. 

[14] For an extremely pedagogical review of the consequences 
of Coulomb frustrated phase separation, see M.Seul and 



D. Andelman, Science 267, 476-483 (1995). 

[15] V.J.Emery, S.A.Kivelson and H-Q.Lin Phys. Rev. Lett. 
64, 475-478 (1990); S.A.Kivelson and V.J.Emery, in 
Strongly Correlated Electronic Materials: The Los 
Alamos Symposium 1993, Eds. K. S. Bedell, Z. Wang, D. 

E. Meltzer, A. V. Balatsky, and E. Abrahams (Addison- 
Wesley, Reading Massachusetts, 1994) pp. 619-656. 

[16] E.L.Nagaev JETP Lett 6, 18 (1967) and 16, 394 (1972); 

P.B.Visscher Phys. Rev. B 10, 943 (1974). 
[17] D.I. Khomskii and K.I. Kugel, Europhysics Letters 55, 

208 (2001). 

[18] There are, in each case, various additional conditions con- 
cerning the nature of the antiferromagnet specified in the 
various conjectural relations between antiferromagnetism 
and Type I behavior. Before gleefully constructing pro- 
posed counterexamples, such as a doped Ising antiferro- 
magnet with additional strong, short-range repulsive in- 
teractions, the reader should consult the original papers. 

[19] It was shown by L. Pryadko, D.Hone, and S.A.Kivelson 
[Phys. Rev. Lett. 80, 5651-5654 (1998)] that a long-range 
uniformly attractive interaction between static charges is 
induced by spin- wave fluctuations in a Heisenberg quan- 
tum antiferrmagnet - these forces are sufficient to prove 
that phase separation is generic in doped antiferromag- 
nets up to a few nasty caveats. 

[20] For a review of stripes in the cuprates see V. J. Emery, 
S. A. Kivelson, J. M. Tranquada Proc. Natl. Acad. Sci. 
USA 96, 8814 (1999). 

[21] C. Howald, P. Fournier, and A. Kapitulnik, Phys. Rev. 
B 64, 100504 (2001); S. H. Pan et al. Nature, 413, 282 
(2001). 

[22] K.M. Lang et al, Nature 415, 412 (2002). 

[23] For a review of the possible role of intrinsic electronic 
inhomogeneity in the manganates see, e.g., E. Dagotto, 
T. Hotta, A. Moreo, Physics Report 344, 1 (2001), and 
Rev. Mod. Phys., in press. 

[24] For a more general experimental review of stripe forma- 
tion in doped antiferromagnets, see J.M. Tranquada, J. 
Phys. Chem. Solids, 59, 2150 (1988). 

[25] B. Keimer et al, Phys. Rev. Lett. 85, 3946 (2000). 

[26] M. Imada, A. Fujimori, Y. Tokura, Rev. Mod. Phys. 70, 
1039 (1998). 

[27] Y. Tokura et al, Phys. Rev. Lett. 70, 2126 (1993). 

[28] Although the cited evidence suggests that LaTiOs is a 
Type II Mott insulator, and thus a counterexample to the 
general conjecture that antifcrromagnetic insulators are 
generically Type I, it seems to us that more work needs to 
be done to establish this fact. The apparent divergence of 
the effective mass with decreasing x reported in [27], and 
the good Fermi liquid character of the metallic state even 
at the smallest reported x, certainly argue in favor of a 



7 



continuous transition to the insulating state. However, 
indirect NMR evidence, Y.Furukawa et al, Phys. Rev. 
B 59, 10550-10558 (1999), apparently implies that there 
is little or no antiferromagnet character to the magnetic 
fluctuations in the metallic phase; a continuous transi- 
tion between a paramagnetic metal and an antiferromag- 
netically ordered insulator would be expected to exhibit 
diverging (quantum critical) antiferromagnetic fluctua- 
tions as the critical point is approached. The failure to 
see such fluctuations is most naturally understood if the 
transition is ultimately first order. The issue is further 
complicated by the fact that the correlation effects are 
relatively weaker in this material than in the cuprates 
or manganates; the charge gap in the insulating state is 
only about O.leV. This might well mean that the range 
of doping over which the actual physics of the "doped an- 
tiferromagnet" is dominant may be rather small, and in 
fact the experiments seem never to access the important 
range of doping with x < 5%. It would be very interest- 
ing to look more closely at the low doping regime of this 
and related materials, and especially to study the mag- 
netic structure by neutron scattering in this regime, to 
unambiguously determine whether they are truly Type 
II. 

[29] M.P.A. Fisher and D-H Lee, Phys. Rev. B 39, 2756 
(1989). 

[30] L. Balents, M.P.A. Fisher and C. Nayak, Int. J. Mod. 
Phys. B 12, 1033 (1998). 

[31] The system considered in this Ref.[4] can be "doped" by 
varying the period of the trap lattice. 

[32] For an introduction to the bose Hubbard model, see, 
for example, R.T.Scalettar et al, Phys. Rev. B 51, 8467 
(1995) and especially, G. G. Batrouni and R. T. Scalet- 
tar Phys. Rev. Lett. 84, 1599-1602 (2000) where issues 
of phase separation are explicitly addressed. 



[33] V. N. Popov, Theor. Math. Phys. 11, 565 (1972). 
[34] D.S. Fisher and PC. Hohenberg, Phys. Rev. B 37, 4936 
(1988). 

[35] M. Blume, V. J. Emery, and R. B. Griffiths, Phys. Rev. 
A 4, 1071 (1971) 

[36] For a discussion of the bicritical point scenario in the 
context of the cuprate superconductors, see S.-C. Zhang 
Science 275, 1089 (1997). 

[37] We note that the vortex superfluid discussed here differs 
from the generic boson superfluid because the net vortic- 
ity is zero. The field theory that describes such "particle- 
hole symmetric" condensate has space and time on equal 
footing (i.e. relativistic). This fact makes the later anal- 
ogy with thermal fluctuating 3D superconductors com- 
pletely appropriate. 

[38] S. Teitel and C. Jayaprakash, Phys. Rev. Lett. 51, 1999 
(1983). 

[39] The field strength F^ v (fj, = 0, 1, 2) of a 2+1-dimensional 
EM field has three independent components. The corre- 
spondence of the these three components with the bo- 
son three-current (p,j) is as follows: Fi 2 <-» p, Fqi <-» 

[40] In the correspondence between a quantum superconduc- 
tor at T = and a thermally fluctuating superconductor 
at finite temperature the magnetic field (F12) and electric 
field (foi,Fo2) in the 2+1 dimensional quantum theory 
are mapped on to B z and (B y , —B x ) respectively. 

[41] M. Tinkham "Introduction to superconductivity", 
McGraw-Hill, New York (1975); P.-G. de Gennes, "Su- 
perconductivity of metals and alloys", Benjamin, New 
York (1966). 

[42] D.R. Nelson, Phys. Rev. Lett. 60, 1973 (1988). 

[43] P.W.Anderson, Science 235, 1196-8 (1987). 

[44] L.Ioffe and A.Larkin, Phys. Rev. B 40, 6941-7 (1989). 



